Spatiotemporal rheochaos in nematic hydrodynamics 
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Motivated by recent observations of rheochaos in sheared wormlike micelles, we study the coupled 
nonlinear partial differential equations for the hydrodynamic velocity and order parameter fields in a 
sheared nematogenic fluid. In a suitable parameter range, we find irregular, dynamic shear-banding 
and establish by decisive numerical tests that the chaos we observe in the model is spatiotemporal 
in nature. 
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Wormlike miceiiar solutions of cetyltrimethylammo- 
mum tosylate (CTAT) show rheochaos for a suit- 

able range of concentrations and temperatures, the stress 
vs. shear-rate curve in 0, 0] shows a plateau signalling 
a rheological phase-coexistence and, in the plateau re- 
gion, the time-series of the shear and normal stresses 
at constant shear-rate show deterministic chaos. This is 
remarkable because it means that the rheological equa- 
tion of state 0, that is, the relation between stress and 
shear-rate, even for macroscopic experimental samples, 
has huge fluctuations about the mean A realistic 
description of the complex dynamics of miceiiar systems 
as well as of other rheological oscillations 0, may re- 



quire variables specific to the system in question 
The studies of 0, 0], however, are significant in 
that they find temporal rheochaos in the dynamics of 
the passively advected alignment tensor alone. They use 
the well-established equations of hydrodynamics for a 
nematic order parameter, outside the regime of stable 
flow alignment, but consider only spatially homogeneous 
states [ljj and hence cannot explore the relation of the 
observed chaos to shear-banding Q. 

In this paper we study numerically the equations of 
the traceless symmetric order parameter for a sheared ne- 
matogenic system, allowing for spatial variation. In the 
parameter range termed "Complex" in the phase diagram 
of the studies |ll], UJl of the spatially homogeneous situa- 
tion, we find spatiotemporal chaos. We give evidence for 
this in Fig. |2I for the dynamic instability of shear bands, 
Fig. |S| for the spatial distribution of "stress drops", Fig. 
^b) for the local phase portrait in the chaotic regime, 
and, decisively, Fig. for the number of positive Lya- 
punov exponents as a function of (sub)system size. We 
now present our model in detail and show how these re- 
sults arise. 

Models of complex rheology generally postulate an 
equation of motion for the stress in addition to that for 
the velocit y fi eld. Our work was motivated in part by 
the failure |l5| to find chaos in the the popular Johnson- 
Segalman (JS) [3 model for the coupled dynamics of the 
stress tensor and the hydrodynamic velocity field. Since 
this model is linear in the stress, and contains bilinear 
products of velocity-gradient and stress, it is natural to 



ask if nonlincarities in the stress are the missing ingre- 
dient. The JS equation is most naturally thought of as 
arising from the underlying dynamics of a local alignment 
tensor, the traceless, symmetric nematic order parameter 
Q, measuring in this case the alignment of the miceiiar 
worms. The stress a OP Il7| d ue to Q is defined via a 
free-energy functional FIofE! H H . Nonlineari- 
ties arise naturally in nematodynamics from terms of or- 
der higher than Q 2 in F. It is these that lead to the rich 
dynamical phase diagram, including a variety of routes 
to chaos, in the work of (1(1 Illlll2j . 

The alignment tensor obeys the equation of motion 
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where I is the unit tensor, u is the hydrodynamic veloc- 
ity field, k = (1/2) [Vu + (Vu) T ] and ft = (1/2) [Vu - 
(Vu) T ] the shear-rate and vorticity tensors respectively, 
r is a bare relaxation time, ao and ot\ are parameters 
related to flow alignment, the subscript ST denotes sym- 
metrisation and trace-removal, and 
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is the molecular field conjugate to Q, for a Landau-de 
Gennes |l8l| free-energy functional F with upto quartic 
terms in Q, with phenomenological coefficients A, B, C, 
Ti, T2, and the simplest bilinears in VQ. In mean- 
field theory, the isotropic-nematic transition occurs at 
A = A, = 2B 2 /9C. As in [H [U E3 , we impose plane 
Couette flow with u = yyx, ignoring modifications of 
this flow by order-parameter stresses, and treating 7 as 
a control parameter, but, unlike in that work, we allow 
spatial inhomogeneity in the y direction. 

As in [l2j . we rescale time by the relaxation time t j A* 
at the isotropic-nematic transition, Q as well by its mag- 
nitude at that transition, and distances by the diffusion 
length made from V\ and t/A*. The ratio T2/T1 of 
the Frank constants is therefore a free parameter which 
we have set to 1 in our study. We choose A — and 
a\ = throughout, and Afe = — (2/\/3)ao and 7 with 
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FIG. 1: Plots showing ai(xo,t) vs ao(xo,t) for tumbling 
(panel T), chaotic (C), onset of aligning (C— +A) and align- 
ing (A) regimes. 

values as sp ecified in the figures. In the orthonormal ba- 
sis 03 T = v^72(zz) S T, Ti = y/t/2{^k - yy), 
T 2 = \/2(xy) ST , T 3 = y/2(xz,) S T, T 4 = \/2(yz)s T , 
we expand Q = ^ afTi and work with the equations 
of motion for the five components a, which follow from 
0. These equations are numerically integrated using 
a 4th order Runge-Kutta scheme with a fixed time step 
(At = 0.001). For all the results quoted here a sym- 
metrized form of the finite difference scheme involving 
nearest neighbors is used to calculate the gradient terms. 
We have checked that our results are not changed if 
smaller values of At are used. We have further checked 
that the results do not change if the grid spacing is 
changed and more neighbors to the left and right of a par- 
ticular site in question are used to calculate the deriva- 
tive. This gives us confidence that the results quoted 
here do reflect the behavior of a continuum theory and 
are not artifacts of the numerical procedure used. We use 
boundary conditions with the director being normal to 
the walls. With this, we discard the 1st 6 x 10 6 timcsteps 
to avoid any possibly transient behavior. We monitor the 
time evolution of the system for the next 5 x 10 6 time 
steps (i.e. t=5000), recording configurations after every 
10 3 steps. We have carried out the study with system 
sizes ranging from 100 to 5000. In the time-series analy- 
sis for the Lyapunov spectrum, we run the simulation till 
t = 20, 000, for a spatial system size of 5000, recording 
data at spatial points at gaps of 10. 

We focus on the parameter region labelled 'C or 'Com- 
plex' in [ill Il2| . where the order-parameter dynamics 
shows temporal chaos. We begin by verifying that the 
inclusion of spatial variation does not alter the phase di- 
agram radically. To this end, we construct local phase 
portraits (plotting different components of the order pa- 
rameter against each other) as a function of the param- 
eters entering into the model, the shear rate 7 and the 
tumbling parameter Afe. Shown in Fig. ^ are the lo- 
cal phase portraits for a particular point xq for various 



shear-rates 7, obtained by holding the value of the tum- 
bling parameter fixed at Afe = 1.25. We have checked that 
the character of the phase portrait (space-filling or oth- 
erwise) remains intact upon going from one space point 
to another though there is no phase coherence between 
two such portraits. A closed curve corresponding to a 
limit cycle is seen in the tumbling region of parameter 
space (denoted by 'T') in the figure, while in the 'C re- 
gion of the phase space it is space filling. When one goes 
away from the 'C region of the phase space to the region 
where the director aligns with the flow the points reduce 
to those on a line and eventually in the aligning regime 
where the director has already aligned with the flow it 
is represented by a point. This assures us that the local 
dynamics in the spatially extended case is similar to that 
of the ODEs of [ll [13 • 

The contribution of the alignment tensor to the devi- 
atoric stress is [H E [U <x OP cx a G - «i(Q • G) ST 
where G, defined in (J2J, is the nematic molecular field, 
and the total deviatoric stress is er OP plus the bare vis- 
cous stress which is a constant within the passive con- 
vection approximation. We therefore look at <x OP alone. 
The spatiotemporal nature of the chaos [2lJ in this sys- 
tem is seen strikingly, if qualitatively, in the space-time 
plots of the shear stress (the xy component of (T OP ), ref- 
ered to as £ henceforth, in Fig. Two other figures are 
instructive, showing the periodic tumbling regime (Fig. 
|3J and a transition from spatiotemporal chaos to sta- 
ble flow alignment (Fig. 0}. The parameter values at 
which these are seen, furthermore, correspond well with 
the phase diagram of [ill Il2j |. Our method therefore 
reveals nicely the rich spatiotemporal structure of the 
phase diagram of this system. In Fig. several distinct 
events, such as the persistence, movement, and abrupt 
disappearance of shear bands can be seen, whereas in 
Fig. El the structure is simply periodic. We find that the 
typical length scale at which banding occurs is a fraction 
of the system size. In order to characterize better the 
disorderly structure of the shear bands in Fig. [3 we look 
at the distribution of band sizes or spatial "stress drops" 
as follows. At a given time (say ti), we define a threshold 
So = 0.8, a little above the global mean (E(y,t)) Vt t, and 
map the spatial configuration to a space-time array of 
±1: S = sgn(S — S ). Fig. EJshows the histogram of the 
spatial length of intervals corresponding to the +state. 
We have considered configurations extending over 2500 
spatial points, and the statistics is summed over configu- 
rations sampled at 5000 times (i.e. i — 1,5000). We see 
that the distribution of band lengths in the spatiotempo- 
rally chaotic regime is fairly broad and roughly exponen- 
tial in shape, and the time-averaged spatial correlation 
function of the stress also decays on a length scale of 
about 5 bands. 

The spatially averaged shear stress shows irregular os- 
cillations in ti me, but analysis of this signal based on the 
techniques of |22j does not give reliable results for esti- 



FIG. 2: Space-time behavior of the shear stress in the chaotic 
regime, 7 = 3.678 and Afc = 1.25. Slice taken from a system of 
size 5000 



mating the correlation dimension v mainly because the 
chaos that we observe is quite high dimensional (embed- 
ding dimension jUI^I m > 10). A very long data train 
is then required for the analysis of the spatially averaged 
time series to yield a correct value of the correlation di- 
mensions. We however establish conclusively that the os- 
cillatory behavior of the stress and the first normal stress 
difference shows signatures of spatiotemporal chaos. In- 
stead of trying to implement the correlation-dimension 
method for our spatially extended problem, we study the 
Lyapunov spectrum |22u24| . Further, instead of studying 
systems of ever-increasing size, we look at subsystems of 
size N s in a given large system of size N, i.e., at space 
points j in an interval iq < j < io + N s — 1 where io is 
an arbitrary reference point. For spatiotemporal chaos 
we expect to find that the number of positive Lyapunov 
exponents grows systematically with N s , which is pre- 
cisely what Fig. El shows. For both figures in Fig. H3 , 
we carry out the procedure with two different reference 
points io and find essentially the same curves. Further- 
more, it has been reported in many studies of spatiotem- 
porally chaotic systems |25j that when calculating the 
subsystem Lyapunov spectrum for increasing subsystem 
size N s , one finds that the Lyapunov exponents of two 
consecutive sizes are interleaved, i.e. the ith Lyapunov 
exponent Ai for the sub-system of size N s lies between the 
ith and (i + l)th Lyapunov exponent of the subsystem of 
size N s + 1. A direct consequence of this property is that 
with increasing subsystem size N s , the largest Lyapunov 
exponent will also increase. This is clearly seen in Fig. 
(b). 

In summary, we have shown that the stress field in 
the equations of sheared, nonlinear nematodynamics, in 
a subrange of the parameter space where flow align- 
ment does not obtain, shows spatiotemporal chaos. Our 
numerical study of this model, which extends that by 
[ToL ITU fl2| for the purely temporal behavior, shows ir- 
regular, dynamic shear banding, and finds that the num- 



FIG. 3: Space-time behavior of the shear stress in the tum- 
bling regime, 7 = 5.0 and = 0.9. Slice taken from a system 
of size 5000. 




FIG. 4: Space-time behavior of the shear stress in the chaotic 
to aligning regime, 7 = 4.1 and Afc = 1.25. Slice taken from a 
system of size 5000. 
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FIG. 5: Spatial distribution of "stress drops" (correspond- 
ing to residence intervals in which the shear stress is above 
a threshold Eo = 0.8) in the Chaotic regime (Fig. 3). In- 
set shows the spatial autocorrelation function averaged over 
t=5000 (circles) with the best fit to the form exp[(- Ay /£)"]. 



4 




FIG. 6: Number of positive Lyapunov exponents (top panel) 
and the largest Lyapunov exponent (bottom panel) as func- 
tions of subsystem size N s , for 7 = 3.678, = 1.25, Em- 
bedding dimension for the time series of each space point is 
10. Diamonds and circles correspond to data taken from two 
different regions of space (io=51 and io=101 respectively, see 
text) . 



ber of positive Lyapunov exponents as well as the largest 
of these grow systematically with system size. Since we 
do not know what parameter values in our model corre- 
spond to the experiments, a detailed comparison to data 
is not possible at this stage. However, it seems likely to 
us that the nonlinear relaxation of the order parameter, 
together with the coupling of nematic order parameter 
to flow, are the key ingredients for rheological chaos in a 
variety of problems. Related issues such as the effect of 
two- and three-dimensional variation, hydrodynamic flow 
beyond passive advection, the routes to spatiotemporal 
chaos and, of course, the role of other degrees of freedom 
such as variable micelle length and lifetime U remain to 
be explored. Meanwhile, we look forward to quantitative 
experimental studies of spatiotemporal rheochaos (2(| . 
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